Spin diffusion in n-type (111) GaAs quantum wells 
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We utilize the kinetic spin Bloch equation approach to investigate the steady-state spin diffusion 
in n-type (111) GaAs quantum wells, where the in-plane components of the Dresselhaus spin-orbit 
coupling term and the Rashba term can be partially canceled by each other. A peak of the spin 
diffusion length due to the cancellation is predicted in the perpendicular electric field dependence. It 
is shown that the spin diffusion length around the peak can be markedly controlled via temperature 
and doping. When the electron gas enters into the degenerate regime, the electron density also 
leads to observable influence on the spin diffusion in the strong cancellation regime. Furthermore, 
we find that the spin diffusion always presents strong anisotropy with respect to the direction of the 
injected spin polarization. The anisotropic spin diffusion depends on whether the electric field is far 
away from or in the strong cancellation regime. 

PACS numbers: 72.25.Dc, 71.70.Ej, 71.10.-w 
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I. INTRODUCTION 

The spintronics with the aim to incorporate the spin 
freedom of the electrons to the traditional electronic de- 
vices to design novel devices has attracted much attention 
in the past decades^ - — For the application of such spin- 
tronic devices, the suitable spin relaxation time and spin 
diffusion/injection length are of essential importance^— 
which requires comprehensive investigation for the thor- 
ough understanding of the spin properties in these sys- 
tems. Actually, the spin relaxation and spin diffu- 
sion/injection have been widely investigated both exper- 
imentally and theoretically for a long period^— 

In n-type III-V zinc-blende semiconductors, the 
D'yakonov-Perel' (DP) mechanism^ which results from 
the spin precession under the momentum-dependent ef- 
fective magnetic field (inhomogeneous broadening^) to- 
gether with any scattering process, is identified as the 
predominant spin relaxation mechanism££ i 15 ' 16 In the 
absence of the external magnetic field, the inhomoge- 
neous broadening is mainly supplied by the spin-orbit 
couplings composed of the Dresselhaus termii and the 
Rashba oneJ£ While the former one is due to the bulk 
inversion asymmetry of the crystal, the latter one origi- 
nates from the structure inversion asymmetry and is tun- 
able, e.g., via the electric field along the growth direction. 
The spin manipulation based on the competition of the 
Dresselhaus and Rashba spin orbit couplings is an in- 
teresting issue in semiconductor spintronics In the 
previous works, it was shown that the electron spin relax- 
ation time in n-type (111) GaAs quantum wells (QWs) 
can be significantly enhanced when the in-plane compo- 
nents of the Dresselhaus term in the vicinity of the Fermi 
surface are strongly canceled by the Rashba termor— 
This suggests the intriguing spin diffusion property in 
this system, which is however still not very clear, to our 
best knowledge. Therefore, the goal of this paper is to 
supply more knowledge for the spin diffusion in (111) 



GaAs QWs. 

Differing from the spin relaxation in time domain, 
the inhomogeneous broadening in spin diffusion is deter- 
mined by Uk = m*[fi(k) + j/jflB]/fc x , 2 - where ft and B 
are the spin-orbit coupling and external magnetic field, 
respectively. According to the previous works, the spin 
diffusion properties are strongly dependent on the de- 
tailed form of the inhomogeneous broadening. Specifi- 
cally, in n-type (001) GaAs QWs where the inhomoge- 
neous broadening is introduced by Dresselhaus spin-orbit 
coupling, Cheng et ai— showed that the spin diffusion is 
enhanced by the scattering but suppressed with the in- 
crease of the temperature in the strong scattering regime. 
In contrast, the spin diffusion length is reduced by in- 
creasing the scattering strength when the scattering is 
rather weak. For the inhomogeneous broadening supplied 
solely by the external magnetic field in Si/SiGe QWs, 
Zhang and Wu 27 found that the spin diffusion is sup- 
pressed by scattering monotonically. Interestingly, the 
spin diffusion length is reduced by the scattering in the 
weak scattering regime and turns to be a constant in the 
strong scattering limit (independent of the electron den- 
sity, temperature, and scattering strength) in the mono- 
layer graphene, where the Rashba spin-orbit coupling is 
dominant^ Due to the partial cancellation between the 
Dresselhaus and Rashba terms, the influence of the elec- 
tron density, temperature, and scattering strength on the 
spin diffusion length in n-type (111) GaAs QWs is ex- 
pected to be an interesting issue. 

In the present paper, we investigate the steady- 
state spin diffusion in n-type (111) GaAs QWs by 
solving the kinetic spin Bloch equations (KSBEs), 5 
which include all the relevant scatterings, i.e., electron- 
impurity, electron-acoustic/longitudinal-optical-phonon, 
and electron-electron scatterings. This work only focuses 
on the DP-limited spin relaxation in the strong scattering 
regime. We find that the electric field (along the growth 
direction) dependence of the spin diffusion length shows 
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a peak due to the strong cancellation between the Dres- 
selhaus and Rashba spin-orbit couplings. We find that 
the spin diffusion length, out of the strong cancellation 
regime, is insensitive to the temperature, doping density 
and electron density. However, the spin diffusion length 
is markedly manipulated by changing temperature and 
doping density around the peak. For the degenerate elec- 
tron gas, the electron density can also lead to observable 
influence on the spin diffusion length. Last but not the 
least, we investigate the anisotropy of the spin diffusion 
with respect to the spin polarization direction of the in- 
jected electrons. The analytical solution of the KSBEs 
with only electron-impurity scattering is also presented 
to explain the numerical calculation. 

This paper is organized as follows. In Sec. II, we in- 
troduce the KSBEs and give an analytical investigation 
for the case with only the electron-impurity scattering. 
In Sec. Ill, the spin diffusion is investigated by solving 
the KSBEs numerically with all the relevant scatterings 
included. A brief summary is given in Sec. IV. 



II. KSBES AND THE ANALYTICAL 
INVESTIGATION 

We start our investigation from n-type (111) GaAs 
QWs under the infinitc-dcpth-squarc-well approxima- 
tion. The well width a is taken to be 7.5 nm and only 
the lowest subband is relevant in our investigation. The 
electrons with their spins polarized along the direction n 
are injected into the QWs at the left boundary (x = 0) 
and diffuse along the x-axis. The right boundary is set to 
be x — L with L much longer than the spin diffusion dis- 
tance. Therefore, the spin polarization of the electrons 
at the right boundary is always set to be zero. Since no 
magnetic field is applied in our model, the spinors pre- 
cess only due to the DP term. Setting the z-axis along 
the growth direction [111], x-axis along [110] and y-axis 
along [112], the total DP term reads 

fi B (k) = 7 (-fc 2 + Mk 2 z ))k y /2V3 - aeE z k y , 
Q v (k) = - 7 (-fc 2 + 4(k 2 z ))k x /2V3 + aeE z k x , (1) 
fi,(k) =y(kl-3k x kl)/s/6. 

Here, 7 and a are the Dresselhaus and Rashba spin-orbit 
coupling coefficients, respectively. E z stands for the elec- 
tric field along z-axis and {k 2 ) = (7r/a) 2 is the average 
of the operator —(d/dz) 2 over the electron state of the 
lowest subband. It is clear that for the particular electric 
field 

E c z ^(Mk 2 z )-k^h/(2V3ae), (2) 

the in-plane components of the Dresselhaus spin-orbit 
coupling are strongly canceled by the Rashba term<22r— 
k 2 here is the average of k 2 over the imbalance of the 
spin-up and -down electrons. 



The KSBEs read^ii 

d t Pk(x,t) = ~ed x ^(x,t)d k!c p k (x,t) - (k x /m*)d x p w (x,t) 
+ d t p k (x, t) | coh + d t p k {x, t) | gcat . (3) 

Here, p k {x,t) are the single-particle density matrices 
of electrons with the in-plane wave-vector k at posi- 
tion x and time t. m* is the effective mass. ^(x,t) 
is the electric potential satisfying the Poisson equation 
Vl^(x,t) = e[N e (x,t) - N Q }/(aK e ) with N e (x,t) = 
Tr[/9k(^, t)} standing for the local electron den- 
sity. Nq is the background positive charge density and 
N e (x,0) = Nq denoting the initial condition. £0 and 
kq are the vacuum and relative static dielectric con- 
stants, respectively. d t pw{x, *)| con = — i[H(k) • cr/2 + 
Ek(a;, t), p-k(x, t)} is the coherent term with Ek(x,t) = 
— J2q Vqpk-q( x i t) being the Hartree-Fock term^S Here, 
a are the Pauli matrices and V q is the Coulomb po- 
tential within the random phase approximation. 30 The 
scattering term dtPk(x,t)\ scg ^ includes the electron- 
impurity, electron-acoustic/longitudinal-optical-phonon, 
and electron-electron scatterings, of which the expres- 
sions can be found in Refs. andl3ll. 

To speculate the properties of the spin diffusion, we 
first simplify the KSBEs with only the elastic impurity 
scattering included in the scattering term. Then, by tak- 
ing the steady-state condition, i.e., dtpk(x,t) = 0, and 
performing the Fourier transformation with respect to 
6>k [with k = fc(cos#k, sm#k)] ; one obtains 

UP^+fc 1 ) - -2p l k /(v k ri)- lk [a z ,p[- 3 +p l + 3 } 
-p k [a + ,p k +l ]+l3 k [cT^p k - 1 ), (4) 

where cr± = (a x ± ia y )/2, v k = f3 k = 

m*[ 7 (k 2 - 4(k 2 z ))/2V3 + aeE z }/2, lk = ijk 2 m* / (2^6), 
and p l k = ^J^de^e-* 16 *. Here, l/r l k = 

m*Ni d9 k [l - cos{Ww)]U 2 /(2n) is the Zth-order 
momentum-relaxation rate. Ni stands for the bulk im- 
purity density and U 2 = f^^-^^^^I^ 
is the impurity scattering potential with |q| = 

y/2k 2 (l — cos6>k). e(q) and the I(iq z ) are the screening 
function 3 ^ and form factor^ respectively. In the strong 
scattering limit, the electron distribution approaches 
isotropy in the momentum space. Therefore, one can 
involve only the lowest two orders of \l\ (= 0, 1) and 
obtain3!i2£ 

d 2 p° k = -2il3 k [a y ,d x p k }+l3 2 [a x , [a x ,p° k }} 

+ fi[<T v ,[<r v ,Pk]]- (5) 

The steady-state spin vector S°.(a;) = Tr[p°(x)cr] then 
can be obtained from Eq. (|5|) together with the boundary 
conditions S°.(0) = S° and S°(oo) = 0. Obviously, S° are 
determined by not only the spin polarization at the left 
boundary but also the electron density and temperature, 
(i) For the injected electrons polarized along the x-axis, 
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i.e., S° = (S , ",0,0) T , the spin polarization vector of the 
steady state is given by 



P(aO = 



(27T) 2 iV e 



ks °k 



dkTr[p k (.T) CT ] 

VI + A 2 sm(\K k \x + 4 
<^ r7 l |(6) 

ci sin(Kkx) 



(ii) For S°, = (0, S^, 0) T , the spin polarization vector is 



P(x) 



/o 27r7V e 



(7) 



(Hi) For Sg = (0,0,S£) T , one obtains 



P(x) 



dk 

2irN e 



C2 sin(_RTfex) 




Vl + A 2 sm(\K k \x- 



.(8) 



Here, A = (8\/2 — ll)/\/7, = arctan(l/A), Z x (&) = 
I,(fc) = (2 v ^+l)/(v / 7+1472| ( 5 fe |), /„(fc) = l/(2|/3 fe |), 

s, and C2 = 



if fe = VTT272/3 fc , ci = - 



(l + V2)\/l+2\/2' 



(20\/2-24)v / i+2\/2 ^ ^y. j QW temperature, these equations 
give the frequency of the spin precession in the space do- 
main, Kk F , and the corresponding spin diffusion lengths 
l x (kF), ly{kp) and / z (fc_p) for the injected spin polarized 
along x-, y- and z-axes, respectively. Here, kF is the 
Fermi wave vector. Since the electrons are in the quasi- 
cquilibrium state in the strong scattering regime, one can 
still treat Eq. (fj| up to the second order and the scatter- 
ing term can be approximately written as p^ 1 /^ 1 for 
the inelastic scatterings. Moreover, from the symmetry 
of the scattering matrices, 1/t° = and \Jt\ = I/t^T 1 
are still tenable. Therefore, the solutions of the spin po- 
larization are still the same as Eqs. (JB])-©. However, the 
carrier exchange effect for the electrons with different k 
due to the inelastic scattering mechanisms, which is not 
included in these equations, results in a unique spatial 
precession frequency.™ As we will show latter, this effect 
can be clearly seen in the vicinity of the cancellation. 



O^ki)' with i = 0, L, satisfy the Fermi distribution for the 
spinors parallel (antiparallel) to the n direction. In the 
calculation, the Rashba coefficient a is taken to be 28 A 2 
fR.ef.l32h. the initial spin polarization at the left boundary 
is set to be 5% and the other parameters are taken from 
Ref.iH. With p k (x,+oo) obtained from the KSBEs, the 
spin polarization signal along the detection direction rh is 
described by P m (x) = J^^ipkix, +oo)a lil \/N e (x, +00), 
which can be well fitted by22 

P m (x) = Aexp(-x/L tih )cos(KtrhX + 4> t )- (11) 

One then defines the spin diffusion length L t m and the 
spatial precession frequency K t m- The acquisition of the 
unique diffusion length and precession frequency from 
Eq. (TlT)) reflects the carrier exchange effect as discussed 
in Sec. II. 

We first investigate the effect of E z by taking both in- 
jection and detection polarization along ^-direction. The 
spatial precession frequency and the corresponding spin 
diffusion length as function of E z are plotted in Fig.[TJ As 
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III. NUMERICAL RESULTS 

To include the electron-electron and electron-phonon 
scatterings, the KSBEs are solved numerically by em- 
ploying the double-side injection boundary conditions 26 
as 



0k(O,t) 



, (for k x > 0), (9) 
(for k x < 0), (10) 



and ^(0, t) — ^(L, t) — since no in-plane static electric 
field is applied. Here, the spins are polarized along n. 



FIG. 1: (Color online) E z dependence of (a) the spatial pre- 
cession frequency and (b) the corresponding spin diffusion 
length at different temperatures and electron densities. The 
impurity density Ni = and the spin polarization is along 
the z-axis. 

shown in Fig.[IJa) , K tz is approximately proportional to 
\E Z — E°\, where the cancellation electric field E c z is cal- 
culated from Eq. © (E c z = 82 kV/cm at 30 K, 80 kV/cm 
at 100 K, and 77 kV/cm at 300 K). This originates from 
the weak momentum dependence of fik out of the strong 
cancellation regime, where k 2 <C 4(/c 2 ) — 2^/ZaeE^/'). In 
Fig 1(b), a peak is clearly seen in the E z dependence of 
the spin diffusion length. One finds that this peak locates 
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just around E°, which suggests that the peak should orig- 
inate from the cancellation effect of the Dresselhaus and 
Rashba terms. Furthermore, it is revealed that the tem- 
perature and the electron density have marginal influence 
on the spin diffusion length for E z far away from E z as 
predicted by Eq. ([8]). However, when E z is around E%, 
f3k is strongly dependent on the moment and Eq. ([SJ fails 
for the inelastic scattering as we explained in Sec. II. In 
this case, one finds that, the spin diffusion length at the 
peak becomes shorter as the temperature increases. The 
underlying physics lies in the fact that, for higher temper- 
ature, the electrons disperse into a broader range in the k 
space, which means that fewer electrons are located in the 
regime with the in-plane Dresselhaus term significantly 
canceled by the Rashba term.— One also notices that the 
spin diffusion length for No = 0.5 x 10 11 cm" 2 (with the 
Fermi temperature Tp « 21 K) is almost the same as that 
for 10 11 cm" 2 (Tp ps 41 K) at room temperature. This 
is because that the inhomogeneous broadening as well 
as the scattering strength (dominated by the electron- 
phonon scattering) is weakly dependent on the electron 
density in the non-degenerate regime at high tempera- 
ture. When the electron gas enters into the degener- 
ate regime, we find that the spin diffusion length can be 
tuned via the electron density. For example, the diffusion 
length changes from 17.5 /im to 15.2 /im as the density 
increases from 10 11 to 4 x 10 11 cm" 2 (Tp 166 K) for 
E z = 79 kV/cm at 100 K. 
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FIG. 2: (Color online) The magnitude of the steady-state 
spin polarization as function of the position at 300 K with 
iV = 10 11 cm" 2 and N a = aNi/N . 

To observe the role of the scattering strength further, 
we introduce the impurities into our GaAs QWs. The 
magnitude of the spin polarization with and without im- 
purities are plotted as function of position in Fig. [2] In 
the absence of the perpendicular electric field, far away 
from the cancellation condition, the additional impurity 
scattering channel only slightly changes the spin diffusion 
length (~ I %) which is consistent with Eq. (|8|). However, 
around E z , the spin diffusion can be significantly influ- 
enced by the high doping density as shown in the figure, 
where L tz changes from 4.0 (3.0) to 3.3 (2.7) /jm with 



E z = 75 (80) kV/cm by introducing the impurities. 

We also investigate the anisotropy of the spin diffu- 
sion by taking the injection and detection polarizations 
along different directions. The position dependence of 
the magnitude of the spin polarization with E z — and 
79 kV/cm are plotted in Fig.[3l where the notation i- 
j means the injected spin polarization along z-direction 
and the detected one along j-direction. Since the signal 
in the x-y, y-x, y-z and z-y configurations are negligible 
as predicted by Eqs. ©-((HJ), we only show the relevant 
components in the figure. For E z = 79 kV/cm, in the 
vicinity of the cancellation electric field, the spin polar- 
izations for the x-z and z-x configurations also vanish due 
to the absence of the spatial spin precession. Interest- 
ingly, it is seen from Fig. [3] that the spin diffusion lengths 
for the x-x and z-z cases for E z — kV/cm are equal 
[L tx = L tz = 0.76 and L ty = 0.51 /im, well consis- 
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FIG. 3: (Color online) The magnitude of the steady-state 
spin polarization as function of the position with different 
polarization directions in the absence of impurity. We take 
T=I00 K and iV = 10 11 cm" 2 . Here, i-j represents the 
measured spin projection along ji-axis while the injected spin 
polarization is along i-axis. 

tent with the relation l x = l z w l.48l y from Eqs. ©-©], 
while the x-x and y-y cases share the same spin diffusion 
length (L tx = Lty = 9.3 /im and L tz = 17.8 /im) for the 
electric field E z — 79 kV/cm. In the latter case without 
spatial spin precession, the drift-diffusion model performs 
well^ Therefore, the spin diffusion length can be ex- 
pressed by L fi = ^DsTsn with r sfl - I/(|f2(k)| 2 - S1?)t p 
standing for the spin relaxation time. D s denotes the 
spin diffusion constant and r p represents the momen- 
tum relaxation timei^ According to Eq. ([T]), one has 
|ft(k)| 2 ~n 2 x = \n(k)\ 2 -tt 2 y « 3.3|ft(k)| 2 -ft 2 , which 
lead to the relation L x — L y 0.55i 2 by considering 
D s and t p free from the small spin polarization. This 
agrees well with the result from the KSBEs (L tx = L ty ~ 
0.52L tz ). 

Finally, we should point out that the correction of the 
envelope function due to the perpendicular electric field is 
neglected, even though it can modify the spin orbit cou- 
pling via the quantity (k z ) and the scattering strength 
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via the form factor in the scattering terni i 30 ! 31 This is 
because: (i) the modification of (k%) due to electric field 
is quite small (within 3 %); {it) the spin diffusion is in- 
sensitive to the scattering strength as discussed above. 
Direct calculation shows that only slight modification of 
the diffusion length is introduced by the electric field, for 
example, around 3 % at 100 K with E z —79 kV/cm. 

IV. CONCLUSION 

In conclusion, we have investigated the spin diffusion 
in n-type (111) GaAs QWs by solving the microscopic 
KSBEs. In the perpendicular electric dependence of the 
spin diffusion length, a peak due to the cancellation be- 
tween the in-plane Dresselhaus spin-orbit coupling and 
the Rashba term is predicted. We find that, for the per- 
pendicular electric field away from the strongest cancel- 
lation value, the spin diffusion length is insensitive to 



the electron density, temperature, and doping density. 
However, in the vicinity of the cancellation electric field, 
the spin diffusion length shows strong dependence on the 
temperature and doping density. For high electron den- 
sity, the electron gas is degenerate, and the spin diffusion 
length is found to be also affected by the electron den- 
sity. Finally, we uncover the anisotropic spin diffusion 
with respect to the spin polarization direction. 
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